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Gas phase mean opacities for varying [M/H], N/O, and C/O 
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ABSTRACT 

We present a set of gas-phase Planck mean and Rosseland mean opacity tables applicable for 
simulations of star and planet formation, stellar evolution, disk modelling at various metallic- 
ities in hydrogen-rich environments. The tables are calculated for gas temperatures between 
lOOOK and lOOOOK and total hydrogen number densities between 10^ cm"^^ and 10'^ cm"-'. 
The carbon-to-oxygen ratio is varied from 0.43 to well above 2.0, the nitrogen-to-oxygen 
ration between 0.14 and 100.0. The tables are calculated for a range of metallicities down 
to [M/H]' = XogNyiINu - -7.0. We demonstrate how the mean opacities and the abun- 
dances of the opacity species vary with C/O, N/O, and [M/H]'. We use the element abun- 
dances from Grevesse, Asplund & Sauval (2007), and we provide additional tables for the 
oxygen-abundance value from Caffau et al.(2008). All tables will be available online under 
http://star-www.st-and.ac.uk/~ch80/datasources.html 
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1 INTRODUCTION 

Star and planet formation modelling (e.g. Wuchterl 2005), stellar 
evolution modelling (e.g. Marques, Monteiro & Fernandes 2008, 
VandenBerg et al. 2008), galactic evolution modelling (e.g. Molla 
et al. 2007) remain dependent on the availability of opacity data. 
Often, the mean opacities were only available on very sparse 
temperature-density grids (Wuchterl 2003; priv. com.). Interme- 
diate mass asymptotic Giant Branch (AGB) stars develop from 
oxygen-rich (C/0<1) into carbon-rich (C/0>1) and are a major 
. source of element enrichment of the insterstellar medium (ISM) 
• (e.g. Gail & Sedlmayr 1986, Dominik et al. 1990, Dorfi & Hofner 
1998, Wachter et al. 2008). However, mass-loss formulae describ- 
ing the matter return from the AGB stars into the ISM do still model 
the gas opacity by a constant values (Wachter et al. 2008).The mod- 
elling of single stars does allow the treatment of a frequency de- 
pendent radiative transfer also in dynamic models for carbon-rich 
(Hofner et al. 2003, Gautschy-Loidl et al. 2004, Nowotny et al. 
2005) and oxygen-rich (Woitke 2006) pulsating AGB stars . Ex- 
tremely metal-poor giant stars with a metallicity of [M/H]< -3.5 
(Cohen et al. 2008, Cayrel et al. 2004) seem to provide strong con- 
strains on the element yields of first super-novae (SN) into the ISM 
(e.g. Cayrel et al. 2004). Therefore, the need for gas-phase opacity 
tables for a wider range of metallicities ([M/H]), carbon-to-oxygen 
ratios (C/O), but also nitrogen-to-oxygen ratios (N/O) has grown 
with the new challenge to model the evolution of star formation at 
high z, i.e. in a comparably young universe with stars containing 
much less elements than our Sun. The actual determination of the 
element abundances profs difficult (see Tsuji 2008) and is particu- 
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lar severe for the solar values (Grevesse, Asplund & Sauval 2007, 
CafFau et al. 2008) which are used as reference throughout astron- 
omy. Furthermore, we need to be aware that stars in other galaxies 
can have an abundance patters different to the Sun e.g. due to a dif- 
ferent star formation history (e.g. Koch et al. 2008). Helling, Win- 
ters & Sedlmayr (2000) presented the difi'erences in the gas-phase 
mean opacities for the small and the large Magellanic cloud which 
element abundances clearly difi'er from that of the Sun. Based on 
this early work on mean opacities for AGB star wind models, we 
provide an extensive grid of mean opacity tables with an increased 
number of line opacity species (Sect.[2j2} for varying [M/H]', C/O, 
and N/O (Sect. [231 and [Sjl. We investigate the influence of these 
parameters on the Planck and the Rosseland mean gas-phase opac- 
ities, and the influence of the oxygen abundance values on our mean 
opacities stimulated by the ongoing discussion (Sect.O. 



2 METHOD 

2.1 Mean opacity formulae 

We compute mean opacity tables in the Rosseland mean approxi- 
mation. 
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with By(T) the Planck function at frequency v and local temper- 
ature T in [K]. «i is the number density in [cm"'] of the gas- 
phase opacity species, p is the local gas density in [gem"']. The 
frequency-dependent absorption coefficient k'^{T) is derived from 
pre-tabulated, opacity-sampled line lists for each of the gas-phase 
opacity source i (i = 1 . . . N). The integrals in Eqs.[T]and|2]are eval- 
uated at 5608 sampling points in frequency space which are dis- 
tributed across the entire wavelength interval 0.2yum . . . 20/jm. This 
distribution samples densest the wavelength ranges were the lo- 
cal Planck functions between lOOOK . . . 6000K peak. These Planck 
functions are meant to cover a range of local temperatures where 
molecules are important opacity sources. For more details see 
Helling & J0rgensen (1998). Helling & J0rgensen (1998) demon- 
strated that the temperature-pressure structure of M-stars is not af- 
fected by the number of sampling points for which the radiative 
transfer is evaluated if changed from ~22 000, to 5608, to ~500 
points. Differences in the temperature-pressure structures become 
apparent for carbon-rich models and for low-metallicity but only 
if the sampling number is as low as ~500. Carbon-rich model at- 
mospheres and low-metallicity model atmospheres show generally 
a higher standard deviation with decreasing number of sampling 
points well below 5000 if looked at the changes in the temperature- 
pressure structure. However, the evaluation of the integrals in Eqs.[T] 
and|2]will to some extend depends on the integral dicretization, and 
hence, on the number of frequency points taken into account. 

2.2 Gas opacity and gas-phase chemistry data 

The number density of the gas-phase opacity species, n\, is com- 
puted for a given temperature, T [K], and total hydrogen number 
density, nn [cm"'], assuming chemical equilibrium for 14 elements 
(H, He, C, N, O, Si, Mg, Al, Fe, S, Na, K, Ti, Ca) and 158 molecules 
with equilibrium constants fitted to the thermodynamical molecu- 
lar data of the electronic version of the JANAF tables (Chase et al. 
1986). The equilibrium constant for TiC are from Gauger et al. (see 
Helling, Winters & Sedlmayr 2000), for CaH from Tsuji (1973), 
and FeH from Burrows (2008, priv. com.). First ionisation states of 
the elements are calculated. The element abundance data are those 
of Grevesse, Asplund & Sauval (2007), and are summarised in Ta- 
ble[T] otherwise will be stated (see Sect. |3.1.2t . The element abun- 
dance are adjusted according to the metallicity [M/H]', the C/O, or 
the N/0 if needed (Sect.[23). 

We include continuum opacity sources (HI (Karzas & Lat- 
ter 1961), H- (John 1988), H+H (Doyle 1968), H^ (Somerville 
1964), H+ (Mihalas 1965), Hel (Peach 1970), He- (Carbon et al. 
1969), and CI, Mgl, AH, and Sil (Peach 1970). Thomson scat- 
tering and Rayleigh scattering for HI and Hel (Dalgamo 1962)). 
The line opacity sources are taken from the Copenhagen SCAN 
data base (CO (Goorvitch 1994), TiO (J0rgensen 1994b), SiO 
(Langhoff & Bauschilder 1993), H2O (J0rgensen & Jensen 1993), 
CH (J0rgensen et al. 1996), CN (J0rgensen & Larsson 1990), C2 
(Querci, Querci & Tsuji 1974), C3 (J0rgensen et al. 1989), HCN 
and C2H2 (J0rgensen 1990 in an updated version), H2-H2 and H2- 
He collision induced absorption from Borysow (2002 and priv. 
com.), and were extended (Schnabel 2001) by the HITRAN data 
base (CH4, NH3, CO2, SO2, NO2, NO, OH, N2; Rothman et al. 
2003). Line list from all these data were opacity sampled according 
to Helfing & J0rgensen (1998). 

In the most general sense, the opacity data can only be as 
good as the line-list data are. To access the quality of line lists, 
studies were undertaken to which extant different line lists can 
help to fit observed spectra. For example, Jones et al. (2002) used 



Table 1. Element abundances, e^i = f^iAsSa _^ from Grevesse, Asplund & 
Sauval (2007) for calculating the gas-phase chemistry. 
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Na 


6.17 


Ti 


4.90 


Mg 


7.53 


Fe 


7.45 



Phoenix model atmosphere structures and different water line lists 
to produce synthetic M-dwarf spectra in comparison to observa- 
tions. J0rgensen (2003) discusses, amongst other molecules, the 
completeness of H20-line list data, and its implications for com- 
parisons as in Jones et al. (2002). In addition to the line lists used 
in this paper, various other line lists are available also for TiO (Plez 
1998, Schwenke 1998) and HCN (Harris et al. 2002, 2008; see also 
J0rgensen 2005). The HITRAN spectroscopic data are regularly 
updated (Rothman et al. 2005) but the data used cover a maximum 
temperature interval of ~ 400 . . . ~ 3000K only. They were ini- 
tially compiled for planetary systems. Much work is done to com- 
pile more complete line lists, one example being CH4 (Borysow et 
al. 2003, Nassar & Bemath 2003, Rothman et al. 2005) for which 
Homeier, Hauschildt Sl Allard (2003) present a line list data com- 
parison. It therefore needs to be appreciated that the line-list date 
we are using are incomplete. We tested the effect of missing weak 
lines and uncertainties thereof (e.g line strength) on the calcula- 
tion of mean opacities. We artificially increased the line opacity 
for lines below a certain threshold. Two cases are shown in Fig.|4l 
(i) all wavelength with k(X) < 10"'cmVg were enhanced by a fac- 
tor of 3, and (ii) all wavelength with k{X) < 10"'cm^/g were en- 
hanced by a factor of 10. These thresholds are small compared 
to the strongest H2O or TiO line of k{A) = 1 - lOcm^/g at e.g. 
Ta; 2000K and p x lO^'g/cm'. Figure|4]demonstrates that possibly 
incomplete line lists have the largest effect on the Rosseland mean 
opacities since it is dominated by transparent spectral regions due 
to the harmonic nature of the averaging process (Eq.[Tll. Therefore, 
uncertainties in weak-line opacities do have a considerable influ- 
ence in particular on the Rosseland mean opacities for a tentative 
line opacity threshold of k{X) < 10"^cmVg. At temperatures below 
lOOOK uncertainties in the line opacities become noticiable already 
at a level of k(X) < 10"^cm^/g. 

2.3 Grid values 

All mean opacities are calculated for a grid of 71 x 71 tempera- 
tures, T, and hydrogen number densities, nn- The values of T run- 
ning from lOOOK to lOOOOK, and ;ih running from 1.0 10"cm"' 
to 1.0 lO^cm"'. We calculate mean opacities for varying input ele- 
ment abundances. We vary the C/O ratio, a value which drastically 
alters the system's chemistry as it passes unity, the N/0 ratio which 
also affects the abundance of the resulting molecular species, and 
the metallicity [M/H]'. Carbon in particular is important here, as it 
may well be from low/intermediate mass stars which pass through 
the AGB phase that the majority of carbon is synthesised (Chiap- 
pini, Romano & Matteucci 2003). of particular interest for mod- 
elling high-mass stars (S. Campbell, priv. com.). 

The abundance of an element is taken to be 6ei = log(A'ei/A'H) + 
12.0 (ie. 6h = 12.0). Then, the ratio of the element in question 
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Figure 1. Oxygen-rich (C/O=0.43) and carbon-rich (C/O=2.0) gas-phase mean opacities for varying nitrogen abundance (N/O=0.14, 5.0, 50.0). Each set of 
four plots corresponds to a different metalUcity, with values [M/H]' = -2.0, -3.0, -4.0, -5.0, -6.0 and -7.0. The data is shown for a hydrogen number density 
HH = 1.931 X 10'5cm-3 (p ~ 10-9gcm-3). 
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Figure 2. Gas-phase mean opacities for the carbon-to-oxygen-rich transition at C/O ~ 1. The left set shows a high metalHcity case ([M/H]' = -2.0), and the 
right a low metalHcity case [M/H]' = -5.0) both for N/O=0.14, 5.0 for ;ih = 1.931 x lO'^ cm"^ 



Table 2. Grid values of C/O, N/0 and [M/H]' for mean opacity tables. 



C/O 0.43, 0.80, 0.90, 0.98, 1.0, 1.1, 2.0, 10.0, 50.0, 100.0 
N/O 0.14, 0.50, 1.0, 2.0, 5.0, 50.0, 100.0 

[M/H]' -2.0, -3.0, -4.0, -5.0, -6.0, -7.0 



(here: C, N or all if metallicity is changed) to oxygen can in turn 
be given as el/0 = 10"^"'"^°. We change the table values (Table 
[U for 6c and to give the desired C/O and N/O ratios, for ex- 
ample setting ec = fo + log(C/0). This also changes the overall 
metallicity. We define the metallicity as [IVI/H]' = logiNu/NH), 
where Nm is the total number density of all the atoms heavier than 
He. Therefore, [IVI/H]' = 0.0 is the metallicity such that the to- 
tal number of metal atoms equals the number of hydrogen atoms. 
Hence, for a desired metallicity each metal must be modified by 
<r = + ^°S(J3) with J3 = 10"^/"!' • (N^'^/NaTK The same 
applies if we vary the C- and the N-abundance (M=C, N) values 
only. Note that our metallicity is not given relative to the solar val- 
ues but is an absolute value. The relative metallicity values, [M/H], 
would write []V[/H] = []V[/H]'-[]V[/H]o with [M/H]o the solar metal- 
licity (see Sect. |3.1.2] ). 

The grid of values for the metallicity and C/O and N/O ratios 
are shown in Table |2] The large number of C/O and N/O values 
should allow models to use the opacity tables for a number of dif- 
ferent purposes, like star and planet foimation or disk simulations. 

We note, that our opacity tables are calculated for a hydrogen- 
rich environment (with the solar abundances as reference) in con- 
trast to the helium-rich environment of e.g. white dwarfs. Hence, 
we have not changed the relative abundance of hydrogen to he- 
lium, and we also refrain from adopting a non-ideal equation of 
state. For example, Harris et al. (2004) demonstrate the influence of 
HeH^ on the opacity of a helium-rich atmosphere. HeH^ is, hence, 
suggested as an important opacity source for metal-free Population 
III stars. Also H^ was appreciated as candidate for efficient cooling 
in metal-free gases which is, however, not included in our mean 
opacity calculation. 



3 RESULTS 

3.1 Opacity tables 

Figure[T]shows plots of several opacity tables for different combina- 
tions of the C/O and N/O ratios, and the metallicity (Table |2j. The 
effects of increasing the N/O ratio and decreasing the metallicity 
are clearly seen to produce an overall reduction in the mean opac- 
ities. The effect due to metallicity can be explained as the fewer 
metals are present, the fewer complex molecules can form, hence, 
decreasing absorption coefficient. The effect of changing metallici- 
ties on the mean opacity values dominates over the changes in C/O 
or N/O. 

One interesting change is seen in the carbon-rich Rosseland 
mean. The two strong opacity peaks in the lower temperature region 
grow smaller as the metallicity is lowered. Eventually, however, 
they, and the trough between them, form a distinctive plateau by 
[IVI/H]' = -5.0. In the ultra-low metallicity cases of [M/H]' - 
-6.0, -7.0, the opacities' features are lessened even more, so that 
by [M/H]' = -7.0 both the Planck and Rosseland means are very 
low and almost constant before the increase due to the continuum 
opacity sources begins to dominate at about log(r5) ~ 3.5. This 
increase at high temperature is seen in all data, but the lower the 
metallicity the lower the temperature at which it can be seen due to 
the drop in molecular contribution to the opacity. 

Interesting to note is that increasing the N/O ratio also ap- 
pears to lower both the Planck and Rosseland mean opacities con- 
siderably. It should however be taken into consideration that by in- 
creasing the nitrogen content, the other metals' abundances will be 
decreased accordingly in order to agree with the given metallicity. 
The changes by increasing the relative nitrogen abundance is com- 
plex because it extensively afi"ects the chemical composition of the 
gas phase (Sect. |T2] l. 



3.1.1 The oxygen-to-carbon-rich transition at C/0~ 1 

The dichotomy between oxygen-rich (C/O < 1) and carbon-rich 
(C/O > 1) opacities is very noticeable. This is mainly due to the 
locking of the majority of either carbon or oxygen, whichever is 
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Figure 3. Differences in Planck mean and Rosseland mean (top panel) opac- 
ities using different oxygen abundance (Grevesse, Asplund & Sauval 2007 
^GrAsSa ^ g gg. c^^a^J^ gt al. 2008 eg'*™ = 8.76). The logarithmic ra- 
tio of the opacities is plotted against temperature. In the lower panel, dif- 
ferences in the mean opacities are also shown for the oxygen abundance 
recommended by Lodders (2003, eo = 8.69), Grevesse & Noel (1993, 
eq = 8.87, used in Helling, Winters & Sedlmayr 2000), and Alexander 
(1975, eo = 8.8093). Note that changing the oxygen abundance results in 
changes of C/O, N/0 and [M/H]'. See text for details. 



the least abundant, in CO. This leaves comparatively little to form 
other molecules, while the more abundant of the two elements may 
do so freely. The changes that occur as C/O passes 1.0 can be seen 
in more detail in Figure|2l Here little difference is seen between the 
C/O ratios of 0.98 (solid lines) and 1.0 (long dashed lines), particu- 
larly in the Planck means. With these ratios, little of either element 
will be available. Meanwhile, at C/O = 1.1 (short dashed lines) a 
large difference is seen as it becomes possible for carbon to form 
opacity-source molecules. Changes in the nitrogen abundances do 
not affect the mean opacities at C/0« 1 . 

Figure |2] demonstrates that the Planck mean values are only 
mildly affected by the drastic chemistry changes in the oxygen-to- 
carbon-rich transition, but the Rosseland mean values do change in 
particular at low temperatures. This is an indication that line opac- 
ity sources with a large number of weak lines (like OH, CH, CN 
or HCN) are important opacity carrier here. The dramatic effect on 
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Figure 4. Comparing mean opacities of this work (red dashes) to Helling, 
Winters & Sedlmayr (2000, blue solid). In comparison, the effect of weak- 
line opacities is shown for two cases (green dots and green long dashes; for 
details see Sect. 12. 2^ 



the atmosphere structure for the oxygen-to-carbon-rich transition 
has been discussed for instance in J0rgensen (1994a). 



3.1.2 The oxygen abundance 

We show how modifying the oxygen abundance affects the result- 
ing gas-phase mean opacities. Two sets of opacity tables were made 
where only the oxygen abundance differ. We test how our mean 
opacities change if we use e^"*" = 8.76 suggested by Calfau et 
al. (2008) in comparison to eO"^A.sSa _ g suggested by Grevesse, 
Asplund & Sauval (2007). Both values are lower than the pre- 
viously used oxygen abundance in Helling, Winters & Sedlmayr 
(2000; see their Table A 1). We use the Grevesse, Asplund & Sauval 

(2007) solar photosphere abundances as specified in Table[T]for all 
other element abundances. Important to note is that while using 
the Grevesse, Asplund & Sauval (2007) value of e^''^'^" = 8.66 
the C/O and N/0 ratios and the metallicity are; C/O = 0.537, 
N/O = 0.132, [M/Hlg'*'^" = -3.02. While using the Caffau et al. 

(2008) values egaffau = § 75^ they are; C/O = 0.427, N/O = 0.105, 
[M/Hlg""™ = -2.97, simply because adding oxygen means adding 
a metal. 

Figure |3] plots for two different gas densities the ratio of the 
mean values for both calculations with different oxygen abun- 
dances. As expected, the mean opacities for e^""™ are higher then 
for e^'^'^'^^'' since more oxygen is available to form molecules which 
are opacity sources. Figure|3](top panel) shows that there is a large 
disparity when changing between the different oxygen abundances. 
While the high-temperature Rosseland and Planck mean values are 
the same since determined by atomic/ionic continuum opacities, 
the deviations can be as large as 20%-30% in the low-temperature 
regime where molecular line absorption dominates the total opac- 
ity. As research into these abundances values is ongoing, tables us- 
ing both values were calculated. 
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Figure 5. Number densities of the line-opacity molecules. The top panels show variation with metallicity, and the bottom variation with N/O ratio («h = 
1.931 X lO''' cm"''). Dashed lines are labelled with their associated values. Solid lines have [M/H] = 3.0 and 4.0 and N/O = 1.0 and 5.0. Each plot contains 
values for C/O=2.0 (red) and C/O=0.43 (blue). Note the scales and ranges of each plot are different. 



3.1.3 Comparing to earlier calculations 

We compare our new mean opacity tables to tables calculated 
with oxygen abundances suggested by earlier works (Fig. [3] lower 
panel): 

fodders ^ g 59 Lodders (2003) 

e'^'^° = 8.87 Grevesse & Noel (1993) 

^Alexander ^ g gQgg Alexander (1975, King mixture IVa), 



and to tables from Helling, Winters & Sedlmayr (2000) for a col- 
lection of Anders & Grevesse (1989) and Grevesse & Noel (1993) 
abundance data (Fig.|4). The oxygen abundance in these old tables 
is 6°"^^" = 8.87 (Grevesse & Noel 1993). For an update of Lodders 
(2003) see Lodders, Palme & Gail (2008). The lower panel in Fig.|3] 
demonstrates that for a given gas density the e^'^" mean opacities 
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Figure 6. Figure |5] is continues liere. Several molecules containing nitrogen are sliown, wliicli demonstrate a dependence on N/O different from the otlier 
molecules. 



deviate by up to 40% in the low-temperature ranges from mean 
opacities calculate using the Grevesse, Asplund & Sauval (2007) 
atmospheric oxygen abundance value. Our mean opacity results 
for the Alexander (1975) eo-value fall in-between the results for 
gLodders e^'^" which is cousisteut with the comparison shown in 
Semenov et al. (2003, their Fig. 3 left). 

The use of the up-dated element abundances by Grevesse, As- 
plund & Sauval (2007) and Caffau et al. (2008), and the use of 



additional opacity sources from the HITRAN data base (Sect. |T2] l, 
causes differences of the present mean opacity values compared to 
our earlier results (Helling, Winters & Sedlmayr 2000). The addi- 
tional HITRAN-data mainly affect the mean opacity values below 
lOOOK. Figure |4] (solid blue vs. red dashed) demonstrates that the 
maximum total difference is 1 - 0.5 orders of magnitude at such 
low temperatures. Semenov et al. (2003, see their Fig. 2) present 
a comparison of our Rosseland and Planck mean gas-opacity data 




Figure 7. Numbers densities for opacity molecules displaying unusual dependence on N/0. Results are only shown for N/0 = 0.14 and 5.0, making it easier 
to see how they vary. HCN and NO2 show that the lines for each N/O may cross one another several times. 



to other works including Alexander (1975). The influence of the 
element abundances on the mean opacities is discussed above with 
emphasis on the uncertainty in the oxygen abundances. Note, how- 
ever, that all element abundances taken from Grevesse, Asplund 
& Sauvale (2008) are lower then in Helling, Winters & Sedlmayr 
(2000, see their Table Al) including the carbon abundances. 

Figure |4] further suggests that possible uncertainties in weak- 
line opacity data can have a considerably larger influence on the 
mean opacities than the uncertainties in the element abundances. 
However, the purpose of our simple weak-line experiment in Fig.|4] 
was to demonstrate the possible influence on the mean opacity cal- 
culation rather then setting a stringent limit on this effects. 

3.2 Chemistry 

In order to evaluate the opacity carriers individual abundance con- 
tribution to the mean opacities, a chemical equilibrium routine was 
run to produce a separate output which gave number densities of 
the various species across the temperature range. The chemistry 
routine considers 158 molecules made of 14 elements, but we only 
present the output for the 18 opacity species included in our mean 
opacity calculation. Varying the C/0 ratio, the N/O ratio, and the 
metallicity, allowed the opacity table's features to be probed (Figs. 
[5]-0. 

The first plot of Fig.|5]shows CO. The density of this molecule 
is effectively constant until log(rg) ~ 3.6, after which it falls away. 
This is true in all cases. Only the height of the plateau and the exact 
point at which it falls off vary with N/O and metallicity. This shows 
why the opacity changes so strongly between oxygen- and carbon- 
rich: Of carbon and oxygen, the least abundant has many of its 
atoms locked within CO, making it harder to form other molecules 
containing that element due to its high binding energy. 

The result of this is seen for many of the molecules shown in 
Figs.|5]and|6l which show large variations between the oxygen- and 
carbon-rich cases. Having C/0 > 1 result in higher densities for the 
carbon-containing species, and the same is true for the case when 
C/O < 1 which affects mostly the oxygen-binding molecules. This 
is entirely as expected. An interesting exception is TiO, where it 
can be seen that at low temperatures (log(r„) ~ 3.1) the densities 
for the carbon-rich case rise above those for the oxygen-rich case, 
because Ti02 becomes more stable than TiO. Lowering the metal- 



licity shifts the point where the densities for each case cross; the 
same occurs when increasing N/O. SiO shows similar behaviour, 
but the rising density as temperature decreases flattens just before 
the high-C density overtakes high-O. Both molecules, TiO and SiO, 
can therefore contribute to the absorption of a carbon-rich gas at 
temperatures below lOOOK. 

Comparing all the plots, similar features become apparent at 
certain temperatures. In the carbon-rich case, it can be seen that 
many molecules' densities have a turning-point at log(rg) ~ 3.4. 
CN, C2 and C3 all reach a peak here before falling again, TiO and 
SiO reach a minimum, and C2H2 and HCN fall off at the end of a 
plateau. The exact temperature at which these changes occur seems 
dependent on the metallicity: in all cases it appears as though de- 
creasing the metallicity shifts the point at which these changes in 
the molecular abundances occur to lower temperatures. 

Looking at the opacity tables for C/0 = 2.0, a peak is no- 
ticed in both the Planck and Rosseland means at log(7'g) ~ 3.4. 
This is caused by the molecules whose number densities also reach 
a maximum here. The peak occurring at very low temperatures 
(log(ri,) ~ 3.1) in the Rosseland mean could be a result of the 
changing densities of several molecules. At the lower temperature 
boundary, where the opacity peak can be seen to be rising, the TiO 
and SiO are at their highest points in the carbon-rich domain, but 
immediately fall. At the same time, C2H2 and HCN are increasing 
to the point where they level out into a plateau. Where they reach 
this point corresponds to the peak in the opacity table. While these 
two molecules stay at roughly the same density, the two oxides are 
removed, lowering the opacity. This holds with Lederer & Aringer 
(2008) who found that C2H2 and HCN are the main contributors to 
their carbon-rich Rosseland mean opacities at these low tempera- 
tures. 

The oxygen-rich molecule densities show a similar 'charac- 
teristic' temperature. At log(rg) ~ 3.4, H2O, which is before this 
point constant, falls off very quickly, while TiO reaches a peak. SiO 
also peaks, but at a point perhaps nearer to log(rg) ~ 3.5. As in the 
carbon-rich case, the metallicity shifts the position of the features. 
The initial plateau in the Planck mean and the prominent peak in 
the Rosseland mean can be easily seen to result from the initial high 
densities of H2O, SiO and SO2, and other molecules which reach 
maxima here such as TiO, OH and NO. 

In the Planck and Rosseland means for both C/0 cases, a small 
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rise can be noticed at logCTg) ~ 3.6. This seems to be due to an in- 
teresting rise in the density of molecules containing the least abun- 
dant of the two elements, eg. if C/O > 1 then molecules contain- 
ing oxygen reach a peak at about this temperature, and vice verse. 
It is important to note that this is the temperature beyond which 
CO's density is seen to fall away, ie. it begins to dissociate. As 
this happens, atoms of both carbon and oxygen are released, al- 
lowing other molecules to form. However the rising temperature 
inevitably brings about continued dissociation of these molecules 
as well, leading to the downwards slope seen on all plots. Some 
of the peaks at log(7'g) ~ 3.6 are of considerable height, allowing 
contribution the rise noted in the opacities. 

3.2.1 Changing nitrogen abundance 

It is also important to note what effects changing the N/O ratio has 
on the molecular densities, and so on the opacities. It can be seen 
that for molecules not containing nitrogen, increasing N/O leads to 
a reduction in the gas density. This is because in order to maintain a 
defined metallicity, when the nitrogen abundance is increased, the 
other metals' abundances must be decreased. For those molecules 
which do contain nitrogen, this is different, and furthermore spe- 
cific to the species in question. For example, the density of NO 
(Figure |6ll increases with N/O in the carbon-rich case, while in the 
oxygen-rich case the opposite is true. Figure|7]shows the difference 
more clearly. Meanwhile, NO2 decreases with N/O in the oxygen- 
rich case, yet when C/O = 2.0 the relationship becomes very un- 
clear with a seeming dependence on temperature as well as N/O. 
Despite these more complex relationships, most molecules still de- 
crease in number with increasing N/O, leading to the overall de- 
crease in the mean opacities observed in Figure [T] 

Comparing the plots of HCN and NO2 in Figures|5]and|6]with 
their counterparts in Figure [T] it may be seen that in the carbon- 
rich cases, and at logCr,) > 3.5 when oxygen-rich, the molecules' 
number densities are lower at N/O = 0. 14 and 50.0 than they are 
when N/O = 1.0 or 5.0. This suggests an optimal N/O ratio which 
is associated with a maximum in the number density of nitrogen- 
containing molecules, and dependent on competition between the 
abundances of carbon, nitrogen and oxygen. Also note that NO2 
remains in very low numbers throughout the whole temperature 
range. While this molecule has little effect on the mean opacities 
in the temperature interval considered here (10' . . . 10*K), HCN 
in particular has a strong influence. 

Finally, also included in the opacity calculations were lines by 
collision induced absorption (CIA) by H2-H2 and H2-He. These 
are known to contribute significantly to lower temperature and high 
pressure opacities, however due to the nature of this kind of absorp- 
tion it cannot be probed as has been done for the molecular lines. 
Borysow, l0rgensen & Zheng (1997) demonstrated that collision 
induced absorption is one of the most important opacity source in 
the atmospheres of the lowest metallicity stars (their Fig. 6). 



4 SUMMARY 

We present a set of gas-phase Planck mean and Rosseland mean 
opacity tables applicable for simulations of star and planet forma- 
tion, stellar evolution, disk modelling at various metallicities of a 
hydrogen-rich gas. The Rosseland mean values are suitable in hot 
and dense regions of an absorbing gas, where the mean free path 
of all photons of arbitrary wavelength is smaller than the charac- 
teristic scale height of the gas, and hence, the radiative transfer 



diffusion approximation holds. Here, the large number of weak ab- 
sorption lines is well represented by the Rosseland mean opacity 
values. If the gas becomes optically thin, the diffusion approxima- 
tion becomes invalid, and single, strong absorbing frequencies may 
become important and determine the radiative flux. In the limiting 
case of a missing correlation between the frequency distribution of 
the flux and the frequency distribution of the opacity, e.g. due to 
line shifts induced by velocity fields, the Planck mean opacity val- 
ues can serves as a good approximation. 

Tables are produced for C/0= 0.43 . . . 100.0, for N/0= 
0.14... 100.0, and [M/H]' = XogNulNn = -2 . . . - 7.0 
([M/H] = [M/H]'-[M/H]o, [M/H]o ^ -3, see Sect l3X2t . We have 
shown that the present uncertainty regarding the oxygen abundance 
values do affect the mean opacity results. We therefor provide ta- 
bles for the values from Grevesse, Asplund & Sauval (2007) and 
for the value from Caffau et al.(2008). We demonstrate the depen- 
dence of our mean opacity values on C/O, N/O and metallicity. The 
oxygen-to-carbon-rich transition range around C/Ow 1 is particu- 
lar interesting, and our results demonstrate that weak line opacity 
sources are important here. We caution, however, that uncertainties 
in the weak-line opacity data can strongly influence the Rosseland 
mean opacities at a tentative threshold of < 10"^cm^/g for k\A) in 
Eq.ID 

The opacity tables will be available online under 
|http://star-www.st-and.ac.uk/~ch80/datasources.html[ 
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